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ABSTRACT 


An  analysis  is  made  of  the  Navy's  demand  forecasting 
process  and  its  impact  on  inventory  system  effectiveness. 

The  current  Navy  Uniform  Inventory  Control  Point  (UICP) 
forecasting  model  is  compared  with  an  alternative  computer- 
oriented  technique  using  UICP  data.  The  comparison  high¬ 
lights  the  presence  of  highly  erratic  patterns  in  the  UICP 
demand  data  base.  Next,  a  simulation  model  is  exercised 
to  suggest  how  the  UICP  demand  reporting  method  might 
contribute  to  the  variance  of  recorded  demand.  The  thesis 
concludes  with  another  simulation  indicating  the  relation 
of  demand  forecasting  accuracy  on  each  component  of  total 
inventory  cost.  This  simulation  suggests  that,  while  holding 
and  ordering  costs  remain  relatively  insensitive  to  fluctua¬ 
tions  in  forecast  accuracy,  the  stockout  cost  element 
displays  a  hypersensitive  reaction. 
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I .  INTRODUCTION 


Inventory  control  is  a  pivotal  activity  of  any  logistics 
organization.  Multi-item  inventory  systems  encompass  trade¬ 
offs  in  balancing  customer  service  needs  with  operating 
costs.  This  management  task  is  particularly  challenging 
in  the  military  setting  where  item  availability  often 
affects  mission  readiness.  Clearly,  one  objective  of  any 
inventory  doctrine  is  to  succeed  in  making  those  decisions 
which  minimize  operating  costs  while  providing  an  acceptable 
level  of  service  for  a  forecasted  rate  of  demand.  The 
demand  forecasting  process  and  its  influence  on  inventory 
system  effectiveness  are  the  subjects  of  this  study. 

First,  some  popular  forecasting  methods  are  profiled  in 
terms  of  six  key  evaluation  measures. 

In  recent  decades,  a  wide  variety  of  forecasting  methods 
has  emerged.  Generally,  they  can  be  assigned  to  one  of  two 
taxonomies.-  qualitative  or  quantitative.  Qualitative 
techniques  are  regarded  as  the  more  subjective  of  forecasting 
approaches.  Usually  conducted  in  a  setting  where  historical 
data  is  unavailable,  this  class  often  employs  expert  opinion 
in  constructing  a  forecast.  The  Navy  Supply  System  uses 
such  a  method,  called  Best  Replacement  Factors,  when  esti¬ 
mating  initial  stock  levels  for  a  new  item  of  inventory. 

In  contrast,  quantitative  methods  make  extensive  use  of 
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historical  data.  There  the  data  serves  as  input  for  various 
types  of  mathematical  models  which  compute  the  required 
forecast,  dot  surprisingly,  advances  in  computer  technology 
tend  to  popularize  the  quantitative-oriented  forecasting 
methods.  Rather  than  spending  considerable  time  synopsizing 
the  more  common  forecasting  techniques,  the  interested 
reader  is  invited  to  consult  the  existing  literature  for 
supporting  detail  (see  Makridakis  and  Wheelwright  [Ref.  1]). 
However,  it  is  interesting  to  classify  a  few  of  the  widely 
used  techniques  in  terms  of  their  cost,  accuracy,  type, 
applicability,  data  pattern  and  time  horizon  characteristics. 
Adapted  from  Wheelwright  and  Makridakis  [Ref.  2],  the  pro¬ 
file  provided  by  Tables  1  and  2  assists  in  selecting  the 
most  appropriate  method  for  a  given  forecasting  requirement. 
For  example,  the  widespread  appeal,  among  industry  and  the 
military,  for  the  exponential  smoothing  technique  is 
apparent.  To  utilize  exponential  smoothing  a  manager  need 
have  only  three  data  elements:  the  most  recent  observation, 
the  most  recent  forecast  and  a  weighting  parameter.  This 
data  storage  consideration  has  been  of  primary  importance 
to  multi-item  inventory  systems  where  demand  forecasts  are 
routinely  prepared  for  several  thousand  items.  Combining 
such  features  as  low  data  processing  and  storage  costs  to¬ 
gether  with  high  applicability,  exponential  smoothing  appears 
as  a  rational  choice  for  the  Navy's  forecasting  method.  In 
contrast,  the  Box-Jenkins  technique  is  not  suitable  for 
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UICP  use.  Although  it  offers  great  accuracy,  the  Box- 
Jenkins  method's  excessive  data  processing  and  data  storage 
requirements  render  its  costs  prohibitive. 

However  profiled  in  terms  of  these  six  factors,  the 
best  measure  of  effectiveness  for  a  forecasting  method  is 
the  economic  benefit  it  provides  to  its  dependent  inventory 
control  system.  The  identification  of  a  method  which  incurs 
low  data  processing  costs  while  contributing  to  minimal  holding 
and  stockout  costs  remains  the  goal  of  the  inventory  fore¬ 
caster.  This  notion  will  be  pursued  further. 

Before  examining  the  cost  effects  of  forecast  accuracy 
several  issues  will  be  analyzed.  The  following  chapter  re¬ 
views  the  forecasting  model  currently  installed  at  one  of 
two  Navy  Inventory  Control  Points,  the  Navy  Ships  Parts 
Control  Center  (SPCC) ,  and  also  introduces  an  alternate 
computer-oriented  model  called,  "Focus  Forecasting."  Chapter 
III  reveals  that  items  having  erratic  demand  patterns  con¬ 
stitute  an  appreciable  portion  of  the  ICP's  inventory 
population.  Chapters  IV  and  V  identify  both  the  cost 
effects  and  some  sources  of  extreme  variability  in  the  UICP 
demands.  Lastly,  Chapter  VI  uses  a  computer  simulation  to 
examine  the  cost  effects  of  forecast  accux-acy. 


II.  TWO  FORECASTING  MODELS 

A.  EXPONENTIAL  SMOOTHING 

As  previously  indicated,  exponential  smoothing  satis¬ 
fies  the  requirements  for  a  forecasting  procedure  which 
demands  low  computer  storage  and  run  time.  Accordingly, 
the  UICP  model  utilizes  a  refined  version  of  exponential 
smoothing  in  establishing,  among  other  things,  the  parameters 
for  the  leadtime  demand  distributions.  Employing  a  modifi¬ 
cation  of  the  continuous-review  inventory  formulation  found 
in  Hadley  and  Whitin  (Ref.  3],  the  UICP  model  applies  the 
forecasted  leadtime  demand  in  determining  both  economic 
order  quantities  and  reorder  points. 

Exponential  smoothing  methods,  originally  advocated  by 
Brown  (Ref.  4],  are  a  geometrically  weighted  sum  of  all  past 
demands  with  the  greatest  weight  applied  to  the  most  recent 
observation.  Mathematically  a  forecast  is  calculated  as: 

NEW  FORECAST  =  LAST  FORECAST  +  WEIGHT  (LAST  OBSERVATION 

-  LAST  FORECAST) 

Notationally  this  may  be  stated  as 


F 


t+1 


Ft  +  a(Dt  -  Ffc)  ;  t  >  1 


which  reduces  to 


t  >  1  . 


F 


t+1 


aDt  +  (l-a)Ft; 


Here  a  is  known  as  the  smoothing  weight  and  is  normally 
assigned  a  value  between  zero  and  one.  The  accuracy  of  the 
exponentially  smoothed  forecast  depends  strongly  on  the 
chosen  a  value.  In  normal  practice,  the  a  smoothing  con¬ 
stant  is  either  selected  arbitrarily  or  suggested  through 
exhaustive  sensitivity  analyses.  With  again  denoting 
the  actual  demand  observation  recorded  in  period  t,  the 
above  equation  can  be  represented  in  recursive  form  as: 


t+1 


=  aD. 


a(l-a)Dt_1  + 


a (1-a)  D 


t-  2 


+  a (1-a)  D 


t-3 


From  this  recursive  form  it  is  easily  seen  how  exponential 
smoothing  dilutes  the  effect  of  the  older  observations. 

The  data-processing  simplicity  of  this  technique  is  evident 
since  one  stored  value,  F.  , ,  replaces  the  entire  block  of 
t  demands.  One  point  of  concern  can  be  the  initial  or  seed 
value  for  Ft .  Since  a  Ffc  value  is  needed  before  the  next 
Ft+^  forecast  is  prepared,  a  seed  value  must  be  found  to 
initialize  the  process  when  t  =  0.  Several  solutions  to  the 
seed  selection  problem  are  available.  Some  of  the  simpler 
solutions  recommend  using  the  first  observation  as  a  seed  or 
dividing  the  time  series  data  into  two  parts  with  the  first 
part  reserved  for  initial  estimation  purposes  (average, 
least  squares  estimation,  etc.).  In  practice  the  seed  problem 
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is  only  of  theoretical  concern.  (See  Makridakis  and 
Wheelwright  [Ref.  1].)  Generally,  the  forecasting  process 
will  be  in  operation  lor.  s  enough  to  suppress  any  dependency 
on  the  seed. 

The  simple  form  of  exponential  smoothing  described  above 
is  essentially  an  approximation  of  a  moving  average  fore¬ 
casting  process.  An  inherent  drawback  of  this  model  is  that 
it  is  relatively  insensitive  to  recent  trend  changes.  This 
weakness,  left  uncorrected,  would  frequently  result  in 
biased  forecasts.  Recognizing  this  limitation,  the  UICP 
model  refines  the  process  by  incorporating  two  types  of 
demand  filters.  (See  Basic  Inventory  Manager’s  Manual  [Ref. 
5].)  First  a  trend  test  is  conducted  to  detect  sustained 
(>_  3)  changes  from  the  underlying  pattern.  The  trend 
statistic,  TR,  is  a  ratio  test  consisting  of: 

2 (SUM  OF  LAST  TWO  OBSERVATIONS) 

TR  (SUM  OF  LAST  FOUR  OBSERVATIONS) 

A  trend  is  considered  present  when  either: 


TR  >  1.1  and  (Dt  ■>  Ft,  Dfc+l  >_  Ffc+1) 

or 

TR  <  0.9  and  (Dt  <  F .  ,  D.  .  <  F  ,  ) 


When  either  of  these  trend  conditions  exist,  the  smoothing 
constant  is  modified  from  its  usual  .10  value  to  a  new, 
"heavier"  weight  of  .30.  This  causes  the  next  forecast  to 
be  more  directly  influenced  by  the  most  recent  observation. 
Next,  a  second  filter  is  used  to  check  for  outlier  observa¬ 
tions.  This  filter  computes  a  control  tolerance  band  around 
using  multiples  of  mean  absolute  deviation  (MAD) ,  the 
average  of  the  absolute  difference  between  actual  and  fore¬ 
casted  demand.  For  further  discussion  on  the  exact  UICP 
use  of  MAD  see  the  Basic  Inventory  Manager's  Manual  [Ref. 

5]  .  If  the  morst  recent  demand  observation  lies  within  a 
band  of  width  7.5  MAD  about  Ffc,  the  process  is  considered  in 
control  and  no  modifications  are  necessary.  If  a  single 
out-of-control  condition  has  been  indicated  by  the  filter, 
the  outlier  demand  is  ignored  and  the  forecast  is  left 
unchanged,  i.e.,  Ffc+^  =  Ffc.  Further,  if  two  successive 
demand  observations  lie  on  the  same  "side"  (high  or  low) 
of  the  tolerance  band,  the  next  forecast  is  calculated  as 
the  average  of  the  two  cohort  outlier  demands,  i.e., 

Ffc+1  =  . 5(Dfc_^  +  Dfc) .  This  condition  is  known  as  a  "step 
increase/decrease. " 

This  completes  the  discussion  of  the  SPCC  exponential 
smoothing  model.  Next,  "Focus  Forecasting"  is  introduced. 

B.  FOCUS  FORECASTING 

Focus  forecasting  is  a  new  forecasting  approach  first 
advocated  by  Smith  [Ref.  6] .  His  concept  requires  the 
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dynamic  simulation  capabilities  of  modern  computers  in 
preparing  each  forecast.  The  methodology  is  somewhat 
inviting  due  to  its  overriding  simplicity.  Focus  fore¬ 
casting  emphasizes  a  straightforward,  flexible  design  in 
acquiring  user  understanding  and  confidence.  The  require¬ 
ment  for  transparent  forecasts  which  are  derived  from  simple 
strategies  and  which  capitalize  on  recent  advances  in  com¬ 
puter  technology  motivated  Smith  in  creating  the  concept. 

The  mechanics  of  the  process  consist  of  four  operations: 
backforecasting ,  selection,  application  and  repetition. 

First,  employing  a  dynamic  evaluation  routine,  the  computer 
identifies  from  a  corpus  of  simple  strategies  the  one  which 
would  have  best  forecasted  the  preceding  period's  demand. 

Next,  this  selected  strategy  is  used  in  preparing  the  up¬ 
coming  period's  demand  forecast.  Lastly,  the  process  is 
repeated  until  "optimal"  strategies  have  been  identified  and 
applied  for  each  inventory  item.  Recent  strides  in  computer 
technology  make  possible  Smith's  procedure  which,  when  imple¬ 
mented  for  large  inventory  systems,  requires  great  internal 
processing  speed.  For  clarity,  a  conceptual  example  of  the 
focus  forecasting  algorithm  is  presented. 

•  Item  Demand  History: 

Time  Period  (Quarters) :12345678 

Units  Demanded  (D^)  :  596  388  527  259  270  363  357  250 

Required:  Ffc  where  t  =  9. 
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Strategy  Corpus: 


STRAT(l) — "LAST  PERIOD  THIS  YEAR" 


Ft  "  Dt-1 


STRAT (2) — "AHEAD/BEHIND  THIS  QUARTER  LAST  YEAR" 


Ft  *  <Dt-lxDt-4'/Dt-5 


Routine: 

Step  1 — Eackforecast  for  t-1,  i.e.,  compute  Ft_^ 
STRAT (1)  =  Fq  =  D? 

=  357  ‘ 

STRAT ( 2 )  =  Fg  =  (D7*D4)/D3 

=  (357  x  259)/527 

«  175 


Step  2 — Select  item — "optimal"  strategy 


DIFF(l) 


DIFF(2) 


ABS ( Dg  -  STRAT (1) ) 
ABS (250  -  357) 

107 

ABS(Dd  -  STRAT  (2) ) 
o 

ABS  (250-175) 

75 


since  75  <  107  item  "optimal"  strategy  is  STRAT (2) 


Step  3 — Apply  selected  strategy 


F9  =  D8XD5/D4 

=  (250  x  270 ) /259 

=  261 


Step  4 — Repeat  for  next  item 
(go  to  step  1) 

Smith  has  implemented  a  seven- strategy  version  for  a  major 
commercial  wholesaler  of  hardware.  There  the  system  pre¬ 
pares  demand  forecasts  for  over  100,000  items  each  month. 
Smith  maintains  that  his  procedure  of  adapting  c  series  of 
forecasting  approaches  to  item  demand  will  significantly 
outperform  a  single-formula  process  such  as  the  UICP  exponen¬ 
tial  smoothing  model.  Before  this  claim  is  examined,  we 
briefly  outline  how  focus  forecasting  might  be  successfully 
applied  to  military  use — specifically  leadtime  demand 
parameterization . 

As  an  initial,  albeit  crude,  attempt  at  implementation, 
six  simple  strategies  derived  from  three  general  forecasting 
categories  are  to  be  used.  The  first  two  strategies  come 
from  the  so-called,  "Naive"  class.  These  are  the  previously 
described  "last  period  this  year"  method  and  its  cohort, 

"this  period  last  year."  (Notationally — STRAT(l)  =  Dt-1? 
STRAT ( 2 )  =  Dt_4.)  The  first  method  acknowledges  trend;  the 
second  seasonality.  The  next  three  strategies  are  selected 
from  the  class  of  moving  averages,  all  of  which  model  hori¬ 
zontal  demand  patterns.  The  moving  averages  will  be  computed 
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as  2,  4  and  8  period  averages.  (Notationally — STRAT ( 3 )  = 

*5(Dt-l  +  Dt-2) '  STRAT(4)  =  ,25(Dt-l  +  •••  +  D  4)  ; 

STRAT (5)  =  .125 (Dt_^  +  ...  +  Dt_g).)  Finally,  as  proposed 

by  Bates  and  Granger  [Ref.  7],  a  composite  forecast  is 

adapted  jn  order  to  improve  forecast  accuracy  by  capturing 

information  from  each  forecast  strategy.  (Notationally — 

STRAT (6)  =  . 2 ( STRAT ( 1 )  +  ...  +  STRAT ( 5 ) ) ) . 

The  next  chapter  compares  the  UICP  exponential  smoothing 

and  focus  forecasting  models  using  empirical  data. 
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III.  MODEL  COMPARISON 


A.  EMPIRICAL  DATA 

To  evaluate  the  two  forecasting  models  under  considera¬ 
tion,  a  nine  year  demand  history  was  retrieved.  The  Opera¬ 
tions  Analysis  Department  at  the  Navy  Fleet  Material  Support 
Office  (FMSO)  followed  a  stratified  sampling  procedure  in 
preparing  a  random  sample  of  522  repairable  and  4530  con¬ 
sumable  items.  Each  of  the  5052  items  was  represented  by 
one  master  data  record  followed  by  several  subrecords.  The 
master  record  contains  descriptive  information  such  as 
national  stock  number,  replacement  price,  etc.  Each  sub¬ 
record  contains  demand  quantity  and  demand  Julian  date 
information.  The  complete  -record  format  can  be  found  in 
[Ref.  8] .  For  convenience,  sequential  dates  ranging  from 
0001  to  3285  replaced  the  Julian  dates  for  the  nine  year 
period.  A  FORTRAN  computer  program  aggregated  the  subrecord 
demand  data  into  36  quarterly  "buckets."  (A  complete  listing 
of  the  principal  FORTRAN  IV  source  codes  used  in  this  thesis 
is  available  in  Appendix  A) .  Negative  demand  quantities 
were  assumed  to  be  the  result  of  customer  cancellations. 

When  negative  demands  were  encountered,  the  quarterly  demand 
balance  was  appropriately  reduced.  When  adjusting  for  such 
cancellations,  however,  the  quarterly  balances  were  forced 
to  maintain  strictly  non-negative  values. 


B.  EVALUATION  CRITERIA 


One  aim  of  a  forecasting  process  is  the  minimization  of 
the  total  forecast  error  incurred  over  time.  When  stochas¬ 
tic  demand  rates  are  involved,  a  natural  assumption  is  that 
superior  inventory  control  requires  very  accurate  forecasting. 
This  assumption  seems  reasonable  since,  in  the  extreme  case 
where  demand  rates  are  deterministic,  existing  inventory 
models  achieve  zero  stockout  cost  and  optimal  ordering  and 
holding  costs.  (See  (Ref.  3].)  Unfortunately,  forecast 
accuracy  lacks  an  absolute  standard  of  measurement.  In¬ 
stead,  the  decision-maker  is  free  to  choose  from  a  wide 
variety  of  evaluation  schemes.  Conceivably,  each  scheme 
could  rank  the  forecasting  models  differently.  This  section 
briefly  describes  the  evaluation  criterion  selected  for 
comparing  the  focus  forecasting  and  SPCC  exponential  smoothing 
models.  Chapter  VI  presents  an  alternate  effectiveness 
measure — the  cost  impact  of  forecast  error. 

Traditionally,  yiared  forecast  errors  (SFE)  have  been 
a  useful  determinant  of  accuracy  with  mean  squared  error 
(MSE)  serving  as  a  popular  choice.  MSE  is  defined  mathe¬ 
matically  as: 


MSE 


1 

n 


n 


l 

t=l 


2 


MSE  is  functionally  related  to  variance  and  also  enjoys 
wide  acceptance  as  a  measure  of  "closeness."  Squaring  a 
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forecast  error  provides  two  advantages.  First,  the  alge¬ 
braic  sign  of  the  error  is  disregarded.  This  avoids  distor¬ 
tion  caused  by  offsetting  positive  and  negative  errors. 

Second,  the  squaring  operation  penalizes  large  forecast 
errors . 

SFE  measures  are  not  without  drawbacks.  For  example, 
what  is  the  preference  ordering  for  two  models  posting 
MSE's  of  81.49  and  75.19?  The  second  forecasting  model 
appears  more  accurate  but  with  what  significance?  Another 
drawback  of  MSE  is  that  it  does  not  facilitate  comparison 
across  different  time  intervals.  These  limitations  not¬ 
withstanding,  three  SFE  measures  will  be  used:  mean,  standard 
deviation  and  the  90t^1  quantile. 

The  next  effectiveness  criterion  is  designed  to  check 
for  biased  forecasts.  Bias  distorts  forecast ‘ accuracy 
through  systematic  overforecasting  or  underforecasting. 

Several  statistics  drawn  from  the  distribution  of  forecast 
errors  help  identify  an  unbiased  forecasting  model.  Any 
reasonably  unbiased  model  will  exhibit  two  attributes. 

First  for  large  samples,  the  forecast  errors  should  have  a 
mean  of  zero.  Second,  as  an  indication  of  symmetry  the 
median  forecast  error  should  approximate  its  mean.  The 
standard  deviation  of  forecast  errors  will  also  be  included 
as  a  measure  of  dispersion. 

A  third  important  consideration  in  determining  forecast 
accuracy  is  how  well  the  model  performs  for  the  higher 
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priced  inventory.  The  manager  may  benefit  little  from  a 
model  providing  low  bias  and  dispersion  for  only  the  lower 
priced  items.  Therefore,  forecast  errors  will  be  weighted 
by  item  replacement  price  (PWFE)  and  tested  by  the  three 
bias  measures  described  above. 

Finally,  the  correlation  between  forecasted  and  actual 
demands  is  considered.  Often  two  time  series  are  closely 
related  but  not  in  a  statistically  dependent  sense.  The 
correlation  coefficient  can  measure  the  strength  of  the 
linear  relationship  between  two  random  variables.  Also 
known  as  Pearson's  product  moment  coefficient  and  denoted 
by  r,  the  correlation  coefficient  is  defined  as: 

n  __ 

l  (D-D)  (F. -F) 

r  =  t=l  C  C _ 

(  l  (D  -D)2  l  (F  -F)  2)  1/2 
t«l  C  t=l  C 

where  D  is  mean  actual  demand;  F  is  mean  forecasted  demand 
r  £  (-1/11. 

A  forecasting  mechanism  which  parallels  changes  in  the 
demand  pattern  is  highly  desirable.  A  high  r  value  would 
typically  indicate  such  behavior  although  it  must  be  re¬ 
garded  with  caution.  For  example,  a  high  r  value  may 
indicate  a  spurious  relationship  caused  by  chance  or  by 
the  elimination  of  a  third  explanatory  variable.  Nonethe¬ 
less,  r  remains  a  recognized  measure  of  the  tracking  capa¬ 
bility  of  a  forecasting  model. 
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In  summary,  four  classes  of  evaluation  criteria  will 
be  reported: 

1)  SFE-- (MSE, SDjPgg)  as  measures  of  closeness, 

2)  FE—  (mean, median, SD)  as  measures  of  bias, 

3)  PWFE — (mean, median, SD)  as  measures  of  price  weighted 

bias , 

4)  CORRE— (r)  as  a  measure  of  association  and  direction. 

Before  reporting  the  results  of  the  model  comparison, 

some  amplifying  remarks  are  offered.  The  SPCC  exponential 
smoothing  model  was  translated  into  the  ES  FORTRAN  source 
code  found  in  Appendix  A.  This  code  contains  all  of  the 
provisions  described  in  Chapter  II.  The  seed  values  were 
computed  as  the  average  of  the  first  four  quarters  of  demand. 
To  dilute  the  effect  of  this  selection,  the  program  was 
stabilized  over  demand  years  two  and  three.  The  demands 
representing  years  four  through  seven  were  examined  in  terms 
of  the  selected  effectiveness  measures.  The  FOCUS  FORTRAN 
source  code  represents  the  previously  described  focus  fore¬ 
casting  model.  For  the  six  simple  forecasting  strategies 
selected,  neither  initialization  nor  stabilization  was  re¬ 
quired.  To  preserve  an  equal  footing,  however,  demand  years 
four  through  seven  were  again  used  for  comparison  to  the 
SPCC  model.  Note  that  each  strategy  must  require  no  more 
than  eight  quarters  of  demand  history.  This  constraint  is 
imposed  since  UICP  data  files  presently  access  only  two 
years  of  demand  history.  Also  note  that  the  FOCUS  program 
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lacks  a  high  demand  filter.  Some  basic  filters  are  offered 
by  Smith  but  none  were  coded  for  the  initial  trial  run. 

The  tables  which  follow  show  the  results  of  testing  the 
SPCC  exponential  smoothing  and  focus  forecasting  models 
using  four  years  of  actual  UICP  demand  history.  To  facili¬ 
tate  model  comparison  between  repairable  and  consumable 
inventory  items,  the  demand  history  was  separated  into  two 
populations.  These  tables  record  each  forecasting  model's 
performance  in  terms  of  the  previously  described  forecasting 
effectiveness  measures. 

C.  FINDINGS 

The  results  in  Tables  3  and  4  do  not  suggest  a  clear 
modelling  preference.  However,  they  do  reveal  some  unexpected 
results  which  led  to  a  redirection  of  this  study.  First 
the  MSE  and  SFE  Pgo  values  indicate  that  the  presence  of 
large  outliers  severely  distorts  MSE.  The  use  of  PgQ  as 
a  more  stable  measure  does  not  indicate  a  significant  dif¬ 
ference  in  model  closeness.  Second,  both  models  succeed 
in  providing  unbiased  forecasts  with  the  focus  forecasting 
model  generally  producing  a  smaller  price  weighted  forecast 
error.  Next,  there  does  not  appear  to  be  any  significant 
difference  in  the  model’s  tracking  ability  as  measured  by 
r.  However,  both  models  appear  more  successful  at  tracking 
consumables  than  repairables . 
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Correlation  r  .98  .98  .98  .99  .98 


The  next  originally  planned  task  was  to  fine  tune  focus 
forecasting  using  high  demand  filters  and  inject  alternate 
forecasting  strategies  for  subsequent  trial  runs.  However, 
this  task  was  interrupted  when  demand  patterns,  such  as  the 
three  actual  demand  histories  shown  below,  were  observed. 


DEMAND  HISTORY 


QUARTER 
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19 
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It  was 

startling  to 

observe  a  significant 

amount 

of 

varia- 

bility  in  both  the  repairable  and  consumable  demand  histories. 
To  gain  a  rough  estimate  of  the  magnitude  of  this  dispersion 
the  average  coefficient  of  variation  (STANDARD  DEVIATION  [ D] / 
EXPECTED  VALUE  [D] )  was  computed.  The  two  data  sets 
recorded  average  coefficients  of: 


Repairable 

COEFF  VAR  = 

3.80 

Consumable 

COEFF  VAR  = 

3.58 

Since  coefficients  of  variation  greater  than  1.0  are  usually 
regarded  as  an  indication  of  excessive  dispersion  [Ref.  1J , 
it  appears  that  the  demand  histories  profiled  above  are  the 
rule  rather  than  the  exception.  Recognition  of  the  data's 
extreme  variability  is  further  suggested  by  GPCC's  1980 
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change  to  their  model’s  high  demand  filter.  The  control 
tolerance  band  was  revised  from  its  previous  F  ±  1.75  MAD 
to  a  widened  range  of  ±  7.50  MAD.  These  erratic  demand 
patterns  challenge  the  foundation  of  those  inventory  models 
(such  as  SPCC's)  which  assume  a  stationary  demand  rate  from 
quarter  to  quarter.  Thus,  two  new  research  directions  were 
suggested:  (1)  identify  the  causes  of  erratic  UICP  demand 

and  (2)  assess  the  cost  impact  of  forecast  inaccuracy  on  the 
type  of  inventory  models  implemented  by  the  ICP’s. 
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IV.  THE  ECONOMIC  BENEFITS  FROM  VARIANCE  REDUCTION 


Before  studying  the  causes  of  erratic  UICP  demands,  it 
is  proper  to  first  identify  the  marginal  benefits  which  can 
be  derived  from  reducing  demand  variance. 

The  stochastic  backorders  model  developed  in  [Ref.  3] 
contains  three  cost  elements:  order,  holding  and  stockout 
costs.  Mathematically  these  are  represented  in  the  total 
annual  variable  cost  equation  as: 

oo 

K(Q,r)  =  AA/Q  +  IC(Q/2  +  r-y)  +  n\/Q[  /  xh(x)dx  -  rH(r)], 

r 

or 


TVC  =  ORDER  +  HOLDING  +  STOCKOUT  COSTS 

For  our  purposes,  K(Q,r)  is  assumed  to  be  both  differen¬ 
tiable  and  jointly  convex  in  Q  and  r.  Under  these  assumptions , 
the  values  Q*  and  r*  which  minimize  K(Q,r)  are  determined 
by  the  methods  described  below.  It  should  be  noted,  however, 
that  joint  convexity  depends  on  the  particular  distribution 
of  h(x).  Brooks  and  Lu  [Ref.  9]  and  Veinott  [Ref.  10] 
show  that  K(Q,r)  can  be  nonconvex  when  h(x)  is  non-decreasing 
and  r  <  y.  (Nonconvexity  may  lead  to  a  failure  in  the 
optimization  technique  which  follows.) 
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Under  joint  convexity  the  optimal  values  Q*,  r*  (where 
0  <  Q*  <  «,  0  <  r*  <  »)  must  satisfy  the  equations: 


3K 

3Q 


3K 

3r 


0  =  li(A  +  irn(r))  +  ~ 

Q 

Q*(r)  =  [ 2 A  (A  +  7rn(rJ)/IC3  1/2 

oo 

C  =  /  (x-r)h(x)dx)  ]  +  IC 

u  r 

=  ^-(-H(r)]  +  IC  (using  Leibniz's  Rule) 

+  H(r*(Q))  =  QIC/ttX 


(note  7T  must  be  sufficiently  large  such  that  H(r*(Q))  <  1)  . 
Computing  a  numerical  solution  for  Q*(r)  and  r*(Q)  poses  a 
problem  of  composite  dependency.  That  is  to  solve  for  Q* 
we  need  to  know  n(r)  and  thus  r.  Secondly,  to  solve  for  r* 
requires  a  knowledge  of  Q.  As  a  practical  and  accurate 
alternative,  a  five-step  numerical  iteration  routine  is 
invoked. 

1)  Check  for  a  unique  solution,  i.e.,  if 

( 2 A  (A+ttu) /IC] 1//2  <  nA/IC  -  unique  Q*,r* 


2)  Assume  n(r)  =  0  and  starting  with  i  =  1, 


1/2 


compute  Q  . 


[2AA/IC] 


3)  Compute  r^  from  H(r^(Q^))  =  Q^IC/ttA  by  consulting 

"tail"  distribution  tables  of  H(x) 

4)  Compute  Q^+1  by  using  r^  to  find  n(r^)  ,  i.e., 

compute  Q^+-j_  =  [ 2 X  (A+nn  (r^) ) /IC] 

5)  Repeat  steps  3)  and  4)  until  Q^+^,r^+^  converge  on  Q^,r^. 

2 

To  determine  the  effect  of  demand  variance  (o  )  on  K 
requires  a  knowledge  of  3K/3a  where 

K(Q,r,a)  -  AA/Q(r,o)  +  IC(Q(r,o)/2  +  r(a)-y)  +  irAn(r,o)/Q(r,a) 

Unfortunately,  developing  a  determinant  form  of  3K/3o  is 
more  complicated  than  the  earlier  task  of  solving  for  Q*  and 
r*.  The  complication  stems  from  the  composite  dependency 
among  a,  Q  and  r.  That  is,  3K/3o  depends  on  3n(r)/3o  which 
depends  on  r  which  in  turn  depends  on  Q.  For  example,  if 
the  leadtime  demand  distribution  is  assumed  normal  with 
mean  y  and  standard  deviation  a  then 

«n|£i  .  |^[  /Vr)»U,c)<3x] 

r 

=  /  “  |f  4><x,o)dx  +  /  (x-r)  — dc 

r  r 

_  2  2 

where  <j>(x,o)  =  (  1//2tto)  exp  [  (x-y)  /2a  ]  and  r  is  the  solution 
to  H(r,o)  =  ICQ  (r  ,o) /ttA  . 

To  remedy  this  complication  and  to  graphically  illustrate 
the  relation  between  K  and  o,  a  sensitivity  analysis  was 
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conducted  using  the  FORTRAN  program  entitled  INVENTOR.  For 
a  set  of  input  parameters  which  remained  fixed  throughout 
the  analysis  1000  realizations  of  a  were  used  in  computing 
the  three  cost  elements  of  K.  The  particular  set  of  input 
parameters  (A  =  500,  X  =  180,  IC  =  15,  tt  =  500,  y  =  90) 
was  selected  to  ensure  joint  convexity  in  Q  and  r.  The 
1000  values  of  a  are  equally  spaced  over  the  interval  [10,350] 
The  upper  and  lower  limits  of  this  interval  were  selected  to 
reflect  the  coefficients  of  variation  found  in  the  UICP 
sample.  The  upper  limit  of  this  range  was  specified  to 
coincide  with  the  average  coefficient  of  variation  (approxi¬ 
mately  3.80)  reported  in  Chapter  III  (i.e.,  350  =  3.80x90). 
Graphical  representations  of  the  sensitivity  analysis  are 
found  in  Figure  1-4.  To  facilitate  comparison  of  the  changes 
in  each  of  the  four  cost  categories,  each  axis  shown  in 
Figures  1-4  was  normalized  in  the  following  manner: 

(1)  Ordering,  holding,  stockout  and  total  variable  costs 
were  each  computed  1000  times  for  the  1000  different 
values  of  o. 

(2)  The  first  value  of  each  cost  (where  o  =  10)  was  used 
as  a  basis  value. 

(3)  Each  of  the  four  sets  of  the  1000  cost  values  was 
then  divided  by  its  respective  basis  to  produce  the 
vertical  coordinates  used  in  the  graphs.  The  1000  o 
values  were  normalized  by  dividing  each  by  o =  10 
thereby  producing  the  horizontal  coordinates  found 
in  the  graphs. 
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These  graphs  highlight,  in  a  general  sense,  the  theoreti¬ 
cal  cost  effects  of  variability  in  leadtime  demand.  Figure 
1  shows  that  as  a  increases  over  the  relevant  range  the 
ordering  costs  developed  in  the  Hadley-Whitin  <Q,r>  model 
decay  in  a  nonlinear  manner.  In  contrast.  Figures  2  and  3 
show  the  corresponding  effects  on  the  holding  cost  and 
stockout  cost  components  of  the  <Q,r>  model.  There  both 
cost  elements  increase  in  a  roughly  linear  fashion  over  the 
relevant  range  with  (3BC/3c)  >  k  *  (3HC/3a)  (where  k  =  4.5) . 
The  reaction  of  total  variable  costs  to  a  has  been  included 
to  show  the  overall  effect.  It  should  be  remembered  that 
the  effects  of  each  of  the  other  three  cost  elements  will 
be  dampened  when  translated  into  total  variable  cost. 

This  sensitivity  analysis  indicates  that  a  study  of  the 
sources  of  leadtime  demand  variance  would  help  reduce  holding 
costs  but  would  also  generate  particularly  beneficial  effects 
on  stockout  costs.  The  identification  of  the  causes  of 
demand  variance  is  the  subject  of  the  following  chapter. 


V.  HYPOTHESIZED  CAUSES  OF  DEMAND  VARIABILITY 


Demand  randomness  results  from  two  sources.  There  is 
the  truly  patternless  variation  that  is  unpredictable  and 
uncontrollable,  but  there  is  also  the  variation  which  re¬ 
sults  from  the  way  data  are  collected  or  reported  during  a 
given  operational  schedule.  This  second  source  of  variation 
is  the  subject  of  analysis  of  this  chapter. 

Erratic  items  constitute  an  appreciable  portion  of  the 
military's  inventory  population.  Silver  [Ref.  11]  addressed 
this  issue  and  it  is  also  demonstrated  by  the  coefficients 
of  variation  found  in  the  UICP  sample  data  discussed  above. 
Silver  maintains  that  imperfections  in  the  wholesale  manage¬ 
ment  information  reporting  system  are  a  primary  cause  of 
erratic  demand.  (FMSO  calls  their  inventory  reporting  network 
the  Transaction  Item  Reporting  (TIR)  System.) 

The  TIR  system  is  a  tri-level  organization  which  sub¬ 
scribes  to  a  policy  of  geographical  responsibility  for  supply 
support.  The  first  level  of  the  TIR's  multi-echelon  network 
is  the  end-user  or  consumer  level  (e.g.,  a  shipboard  supply 
department) .  Level  two  is  composed  of  several  retail  out¬ 
lets  such  as  the  at-sea  Mobile  Logistics  Support  Force 
(MLSF)  ships  and  the  shore-based  stock  points  (Naval  Supply 
Centers,  industrial  rework  facilities,  Naval  Air  Stations, 
etc.).  The  final  tier  contains  the  Navy's  two  wholesale 
activities-- the  Inventory  Control  Points. 
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Under  the  most  common  scenario,  a  shipboard  supply 
department  receives  a  routine  demand  request  from  one  of 
its  customer  departments.  If  the  demand  reduces  the  ship's 
available  stock  past  the  reorder  point  or  if  the  demand 
is  for  an  unsupported  item,  a  replenishment  requisition  is 
submitted  to  the  nearest  retail  outlet  by  one  of  two  means. 

For  ships  operating  at  sea  with  an  assigned  MLSF  unit, 
the  replenishment  requisition  is  submitted  directly  to  the 
MLSF  ship.  For  ships  operating  without  MLSF  support  or  for 
those  which  are  inport  when  the  demand  occurs,  the  replenish¬ 
ment  requisition  is  deposited  with  the  nearest  shore-based 
retail  outlet.  Clearly  for  situations  where  a  customer-ship 
operates  at  sea  for  extended  periods  without  MLSF  support, 
the  replenishment  requisition  encounters  substantial  delays 
before  reaching  the  TIR's  retail  level.  Under  either  case 
when  the  requisition  ultimately  reaches  a  retail  activity, 
the  TIR  system  activates  and  provides  transaction  status  to 
all  levels.  When  the  available  stock  at  the  retail  outlet 
reaches  its  reorder  point,  a  collective  replenishment  requi¬ 
sition  is  submitted  to  the  ICP .  This  requisition  reflects 
the  current  and  anticipated  future  demands  from  the  retail 
outlet's  geographically  assigned  customers.  In  this  cas¬ 
cading  fashion  every  demand  which  originates  at  the  consumer 
level  ultimately  reaches  the  wholesale  level  where  the  ICP 
i3  tasked  with  the  system-wide  reorder  and  budgetary 
responsibilities . 
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A  simulation  was  conducted  to  illustrate  how  the  TIR 
organizational  structure  may  cause  truly  nun-erratic  con¬ 
sumer  level  demands  to  be  observed  as  the  highly  erratic 
demand  actually  seen  at  the  ICP  level.  (See  Lewis  and  Uribe 
[Ref.  12]  and  FORTRAN  source  code,  DEMSIM.)  The  framework 
for  this  simulation  was  adopted  from  Shields  [Ref.  13]. 

The  simulation  modified  the  present  TIR  demand  reporting 
scheme  by  inserting  an  artificial  reporting  channel  which 
reports  all  consumer  demands  directly  to  the  ICP  as  they 
occur.  To  date  this  reporting  medium  has  been  reserved  for 
only  high  priority  requisitions.  This  alternate  reporting 
channel  provides  two  valuable  capabilities.  First  by  directly 
reporting  user-demands,  the  ICP*s  records  will  accurately 
reflect  the  demand  patterns  occurring  on  the  consumer  level. 
Second,  the  ICP  recording  operation  is  expedited  by  elimin¬ 
ating  the  reporting  delays  incurred  when  ships  operate  with¬ 
out  MLSF  services.  Together  these  accuracy  and  timeliness 
considerations  can  have  a  big  impact  on  the  effectiveness 
of  an  ICP. 

The  structure  for  the  DEMSIM  program  is  as  follows: 

1)  assume  a  30  ship  population  with  each  ship  operating 
independently  over  an  13  year  active  service  life. 

2)  initial  ship  inport  and  at-sea  assignments  were 
randomly  selected.  (The  terms  "inport"  and  "at-sea" 
are  meant  only  to  describe  when  a  retail  supply 
activity  is  readily  available) . 
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3)  assume  inport  periods  follow  a  Uniform  [5,30]  dis¬ 
tribution  and  at-sea  periods  follow  a  Uniform  [5,45] 
distribution. 

4)  generate  demands  in  a  random  manner  following  a 
compound  Poisson  process  where  demand  interarrival 
times  are  Exponential  and  demand  requisition  sizes 
are  Poisson. 

5)  record  demands  according  to: 

a)  the  indirect  reporting  system  presently  used 
by  ICP. 

b)  the  artificial  system  where  the  ICP  receives 
observations  directly  from  the  consumer. 

6)  aggregate  each  demand  record  into  quarterly  groupings 
and  compute  the  variance  of  each  time  series . 

7)  repeat  the  process  50  times. 

The  result  of  the  simulation  supports  a  similar  finding 
to  that  of  Shields  [Ref.  13],  that  demand  variance  is  sub¬ 
stantially  reduced  using  the  alternate  reporting  method. 
After  50  replications  the  direct  reporting  scheme  demon¬ 
strated  an  average  variance  of  3050  as  opposed  to  5560  for 
the  indirect  ICP  reporting  system  la  45%  reduction] . 

This  simulation  is  meant  to  highlight  only  one  cause  of 
erratic  demand.  Presumably  other  changes  in  the  demand 
recording  system  could  yield  additional  improvements. 
Examples  of  such  changes  are: 

*  Limited  Demand  History — continued  strides  in  computer 
technology  will  make  storage  of  longer  demand  histories 
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an  economic  reality.  Greater  depth  in  the  demand 
history  would  facilitate  aggregation  by  other  than  a 
quarterly  basis.  A  change  in  the  base  time  period 
could  provide  a  smoother  stream  of  observations. 

*  Customer  Identity — The  interaction  of  large  customers 
(e.g.,  naval  bases)  who  order  more  product  less  often 
with  small  customers  (e.g.,  destroyers)  who  order 
less  product  more  often  might  cause  high  variability. 

It  may  prove  advisable  to  monitor  a  few  large 
customers  on  an  individual  basis. 

*  Customer  Buying  Habits — A  practice  of  batch  ordering 
could  easily  be  a  contributing  factor.  Actual  consump¬ 
tion  of  material  may  follow  a  stable  pattern,  however, 
due  to  3uch  factors  as  funding  and  increased  operational 
commitments  material  is  often  ordered  in  large  lots. 

*  Product  Identity — The  disaggregation  of  products 
into  their  life  cycle  states  (growth,  steady  state, 
gradual  phase-out)  might  produce  different  and  dis¬ 
tinct  demand  patterns. 

Whether  such  changes  could  or  even  should  be  implemented 
largely  depends  on  the  results  gleaned  from  comparing  the 
inventory  cost  savings  against  the  cost  of  altering  the 
demand  recording  system.  This  study  does  not  undertake 
such  a  task  but  focuses  instead  on  a  more  pertinent  issue — 
the  dollar  cost  of  forecast  accuracy.  That  is,  for  a  given 
demand  pattern  how  desirable  is  it  from  an  inventory  cost 
perspective  to  pursue  a  more  accurate  forecasting  model? 
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VI.  COST  IMPACT  OF  FORECAST  ERROR 


The  final  phase  of  the  study  uses  a  computer  simulation 
to  examine  how  forecast  error  affects  ordering,  holding  and 
stockout  costs.  This  examination  does  not  suggest  how  to 
improve  forecasts  but  it  does  indicate  if  it  is  worth  doing 
anything  at  all. 

The  <Q,r>  model  presented  in  Chapter  IV  was  a  repre¬ 
sentation  of  a  theoretical  inventory  cycle  for  which  sta¬ 
tionary  annual  demand  rates  were  applicable.  Here  the 
interest  is  in  measuring  the  actual  inventory  costs  which 
an  ICP  incurs  when  annual  demand,  due  to  forecasting 
inaccuracy,  is  no  longer  stationary  and  when  both  the 
reorder  quantities  (Qi)  and  reorder  points  (r^  vary  for 
each  ith  cycle.  To  conduct  such  an  analysis  requires  a 
type  of  "bookkeeping"  computer  routine  capable  of  both 
maintaining  the  three  inventory  cost  "accounts"  and  com¬ 
puting  and  r^  for  each  cycle. 

For  the  purposes  of  this  chapter,  the  i  cycle  length 

(T . )  is  defined  as  the  interval  of  time  between  receipt  of 
i 

successive  orders.  Registering  cycle  order  costs  reduces 
to  the  simple  task  of  posting  the  cost  per  order  (A)  once 
each  cycle.  The  second  account,  holding  costs,  is  posted 
in  one  of  two  ways.  First,  if  no  backorders  occur  then  the 
cycle  holding  costs  are  jT^  (Opening  Inventory)  i  +  (Closing 
Inventory).^  ^]  .  If  a  backorder  does  occur,  however,  the 


cycle  holding  costs  are  just  ^-[T|  (Opening  Inventory)^]  where 


Tj  is  that  increment  of  for  which  the  inventory  balance 
was  positive.  It  should  be  noted  that  in  exercising  this 
actual  inventory  cycle  the  opening  cycle  inventory  may 
differ  from  the  Q+(r-y)  quantity  found  in  the  theoretical  inven 
tory  cycle  of  Chapter  IV.  A  more  accurate  representation  is 
nQi+ ( r i”h i )  where  n  is  an  integer.  This  is  necessary  to  ensure 
that  the  reorder  point  for  each  successive  cycle  will  be 
reached.  That  is,  since  and  r^  change  with  each  cycle 
under  certain  situations  an  integer  multiple  of  must  be 
ordered  during  the  i^1  cycle  to  enable  the  i+1  opening 
cycle  inventory  to  exceed  its  reorder  point,  ri+^. 

As  to  the  third  account,  stockout  costs,  two  possibili¬ 
ties  must  be  considered.  The  first  possibility  is  that  the 
backorder  cost  is  rr  per  actual  backorder  with  no  penalty 
assigned  for  the  length  of  time  the  backorder  exists. 

A 

Second,  the  actual  backorder  cost,  tt  ,  could  also  be  con¬ 
sidered  as  a  function  of  the  time  the  backorder  exists 
(Ti  -  T  j) .  Separate  accounts  are  used  to  examine  the  effects 
of  forecast  error  on  each  possibility. 

Before  presenting  a  summary  of  the  structure  of  the 
computer  simulation,  a  measure  of  forecast  error  is  required 
to  demonstrate  the  cost  effects  for  various  magnitudes  cf 
forecast  error.  The  difficulties  of  dealing  with  a  single 
criterion  of  forecast  accuracy  was  addressed  earlier.  In 
spite  of  the  fact  that  it  can  at  times  be  an  unreliable 


6 


figure  of  merit  (see  Tables  3  and  4) ,  MSE  was  selected  as 
the  measure  of  accuracy.  Specifically,  an  array  of  eleven 
equally  spaced  values  ranging  from  0  to  5000  was  chosen 
(i.e.#  MSE ( k)  e  [0,5000],  k  =  1,11).  This  range  of  values 
partially  reflects  the  MSE  quantities  found  in  Tables  3 
and  4 . 

The  simulation  computer  program  (see  MSESIM  FORTRAN) 
borrowed  the  input  parameters  found  in  the  INVENTOR  source 
code  discussed  in  Chapter  IV  (with  a  =  10  here) .  In  the 
MSESIM  routine  a  demand  value  (\^  -  NORMAL (180 ,225) )  and  a 
leadtime  (LT^  ~  NORMAL ( 6 (mos) ,1) )  were  generated  for  each 
cycle  using  an  available  random  number  generator  (see  [Ref. 
12]).  Using  the  synthetic  demand  rate,  the  available  cycle 
stock  level  is  depleted  until  the  cycle  reorder  point  (r^ 
is  reached.  This  triggers  the  placement  of  an  order  of  size 
nQ^  such  that  this  replenishment  quantity  will  exceed  r^+1» 
Stocks  continue  to  deplete  at  a  rate  of  over  the  syn¬ 
thetic  leadtime  period  (LT^) .  Upon  reaching  the  end  of  LT^ 
the  cycle  T^  is  completed  and  the  reorder  quantity  nQ^ 
arrives.  At  this  point  the  accounts  for  ordering,  holding 
and  stockout  costs  are  updated  to  reflect  the  costs  actually 
incurred  during  the  cycle.  Next,  and  r^+^  are  com¬ 

puted  using  the  basic  algorithm  found  in  Chapter  IV  except 
that  A^  is  deliberately  falsified  to  reflect  the  forecast 
error  introduced  by  the  chosen  forecasting  model.  (Fore¬ 
cast  errors  (FE^)  were  randomly  generated  as  NORMAL (MSE (k) , 
1)  deviates.)  In  this  way,  the  FE^  provided  by  the 
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forecasting  model  causes  a  nan-optimal  and  r^+^  to 

be  computed.  This  process  continues  until  a  time  counter 
reaches  a  specified  upper  limit  at  which  time  the  evolution 
is  rerun  for  another  pass.  After  fifty  passes  the  next 
candidate  forecast  error  (MSE(K+1))  is  introduced  and  the 
entire  simulation  process  repeats  for  another  set  of  fifty 
passes.  After  all  eleven  values  of  MSE(K)  are  used  the 
computer  program  plots  the  grand  average  for  each  cost 
account  as  a  function  of  its  corresponding  value  of  MSE(K) 
producing  the  graphs  which  follow. 

The  figures  presented  below  are  the  results  of  the  MSESIM 
computer  routine  just  described.  Again  each  ordinate  axis 
is  normalized  to  facilitate  comparison.  Also  for  complete¬ 
ness  both  forms  of  stockout  costs  (unit  and  time  weighted) 
and  their  corresponding  total  variable  cost  graphs  follow 
the  ordering  and  holding  cost  plots. 

A  review  of  Figures  5-10  shows  that  both  ordering  and 
holding  costs  exhibit  relatively  minor  changes  over  the 
relevant  range  (9%  decrease;  6%  increase,  respectively) . 

The  roughly  linear  increases  in  both  stockout  costs  are 
quite  dramatic  especially  for  time-weighted  backorders  where 
greater  than  a  five  fold  increase  is  observed.  Figures  8 
and  10,  the  two  total  variable  cost  graphs,  must  be  viewed 
with  caution  for,  once  again,  their  behavior  is  directly 
influenced  by  the  particular  values  selected  for  A,  IC,  it. 


Figure  8.  Total  Cost  (Unit  Weighted)  Sensitivity  to  MSE 
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The  conclusions  gleaned  from  these  graphical  displays 
are  two.  First,  the  consequences  of  even  severe  forecast 
error  on  both  ordering  and  holding  costs  appear  relatively 
insignificant.  The  second  and  more  profound  finding  is 
that  a  particular  forecasting  model’s  accuracy  (as  measured 
by  MSE)  appears  to  play  a  vital  role  in  keeping  total 
stockout  costs  reasonable. 
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VII.  CONCLUSIONS 


A.  SUMMARY  AND  FINDINGS 

The  study  began  with  a  comparison  of  the  SPCC  forecasting 
process  to  an  unrefined  form  of  a  newer  computer-oriented 
technique.  The  use  of  several  effectiveness  measures  showed 
neither  exponential  smoothing  nor  focus  forecasting  per¬ 
formed  significantly  better  in  forecasting  actual  SPCC 
demands.  One  outgrowth  of  the  comparison  was  the  prepon¬ 
derance  of  erratic-demand  patterns  in  the  IC?  inventory 
population.  Next  a  sensitivity  analysis  was  designed  to 
measure  the  cost  effects  of  demand  variability.  This  analy¬ 
sis  indicated  that  reductions  in  demand  variability  gener¬ 
ate  generous  cost  savings,  particularly  in  the  form  of 
reduced  stockout  costs.  Not  content  with  viewing  the  erratic 
demand  situation  as  inflexible,  two  research  directions  were 
pursued:  (1)  a  change  to  the  demand  pattern;  (2)  the  cost 

effects  of  changing  the  forecasting  model.  First  the  in¬ 
clusion  of  an  artificial  demand  reporting  channel  produces 
considerable  reductions  in  demand  variance.  A  second  and, 
presumably,  more  complicated  improvement  to  the  demand 
forecasting  process  would  be  the  development  of  a  model 
capable  of  accurately  forecasting  erratic  items.  To  gain 
insight  into  the  expected  economic  worth  of  any  such  model, 
the  study  concluded  with  an  examination  of  the  cost  effects 


of  forecast  error.  There  it  was  found  that  stockout  costs 
decrease  drastically  as  forecast  accuracy  improves. 

B.  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Follow-on  research  in  the  areas  listed  below  may  aug¬ 
ment  the  research  conducted  in  this  study. 

1)  Conduct  further  investigation  into  modifying  the  TIR 
reporting  system  by  routing  all  consumer  demands 
directly  to  the  ICP . 

2)  Investigate  other  causes  of  demand  variability  and 
examine  the  corrective  procedures  necessary  to 
smooth  input  demand  data.  (Several  possible  sources 
were  hypothesized  in  Chapter  V.) 

3)  Continue  to  analyze  computer-oriented  forecasting 
techniques.  Since  the  development  of  the  exponential 
smoothing  concept,  technological  advancements  in  the 
computer  field  have  lessened  computational  and  data 
storage  costs.  Newer,  innovative  forecasting  methods 
capitalize  on  these  advancements.  Research  into 
dependable  forecasting  models  which  are  easily  under¬ 
stood,  easily  maintained  and  easily  adapted  by  a  multi¬ 
item  inventory  activity  should  receive  continued 
sponsorship.  The  focus  forecasting  concept  is  one 
possible  candidate. 
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COMPUTER  PROGRAM  SOURCE  CODE 
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